Quantum Control via Geometry: An explicit example 
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We explicitly compute the optimal cost for a class of example problems in geometric quantum 
control. These problems are defined by a Cartan decomposition of su(2") into orthogonal subspaces 
[ and p such that [[, [] C p, [p, I] = p,[p,p] C [. Motion in the [ direction are assumed to have 
negligible cost, where motion in the p direction do not. In the special case of two qubits, our results 
correspond to the minimal interaction cost of a given unitary. 
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I. INTRODUCTION 

Characterizing the difficulty of synthesizing particular 
quantum interactions has generated considerable inter- 
est in recent years due to its practical applications in 
quantum computation. From the perspective of optimal 
control, it determines the optimal way to construct a 
desired quantum interaction with a limited set of tools 
P, [3, H, |j, Q . From the perspective of quantum circuits, 
it expresses the minimal number of basic gates required 
to build up a given algorithm Q . 

These perspectives result in different characterizations 
of complexity. In optimal control, a unitary is hard if 
it is costly to synthesize with available interactions. In 
quantum circuits, a hard unitary requires a large number 
of the basic available gates. Recent work by Nielsen et 
al. shows that for certain control problems, both charac- 
terizations are polynomially equivalent^ 121]. 

This equivalence motivates the application of contin- 
uous geometrical methods to quantum circuits, and in 
the special case where the basic gates are single and two 
qubit unitaries, quantum complexity [1, Q. In this for- 
mulation, each unitary operator corresponds to a point in 
a particular Riemmannian manifold. The metric is engi- 
neered such that the minimal distance between a unitary 
operator U and the identity / corresponds to the minimal 
cost of synthesizing U. This approach allows us to apply 
mathematical techniques cultivated over many decades 
to a significantly newer field. 

Prior work in quantum optimal control has mostly 
dealt with systems that evolve under a specific drift 
Hamiltonian (See for example, [O, [3 ) ■ However, 
all entangling operations are equivalent modulo local in- 
teractions [13, H, [IE [3 [i3| and hence no particular op- 
eration should be favored in a model compliant with the 
spirit of quantum complexity. This motivates the treat- 
ment of interaction Hamiltonians as a physical resource, 
where they are all assigned equal cost. 



In this paper, we consider a class of quantum control 
problems where the space of Hamiltonians is divided into 
two orthogonal subspaces, the application of Hamiltoni- 
ans in one subspace incurs negligible cost compared to 
the other. Provided these subspaces satisfy the condi- 
tions of a Cartan decomposition (see below), geometrical 
methods may be used to construct a general solution. 

While the general class of systems solved in this pa- 
per has not been analyzed in previous literature, it en- 
compasses a number of previously studied systems. In 
the special case of a single qubit, our result provides an 
alternative characterization of single qubit time-optimal 
control In the case of 2-qubits, our solution coin- 

cides with the interaction cost of a two qubit unitary , 
minimized over all possible drift Hamiltonians. 



II. BACKGROUND AND DEFINITIONS 

In this section, we introduce some of the necessary 
background and notation that will be used in the paper. 
We assume the reader is familiar with the basic notions of 
Riemannian geometry. Lie algebras and quantum circuits 
(E.g 6, 19, 20]). 

Consider an n-qubit system. The space of traceless 
Hamiltonians H S su{2") on this system forms a vector 
space under the trace inner product {A\B) = tr {AB). 
This space is spanned by the product operator basis 
rijLi where crj,k denotes the action of applying the 
Pauli interaction cTj e {/, ax, cry, cFz] to the fc*'' qubit. 

The quantum control problem is defined as follows. We 
wish to synthesize a given n-qubit unitary U € SlUT^) 
by the application of some Hamiltonian H{t) S sm(2"). 
We define a cost function C : S'C7(2") x su(2") M such 
that the application of H for duration dt on a unitary U 
incurs cost C{U,H)dt. Formally, the system is governed 
by the Schrodinger equation 



'dt 



-iH{t)U(t) U{0) = / U(T) = U. (1) 



^ The are some technical caveats to this equivalence related to 
approximate versus exact implementation. See [3 for details 



We aim to find the H{t) on interval [0, T] such that the 
total cost D{I, U) = Jp C{U, H)dt is minimized. 
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In this paper, we analyze such problems using the ge- 
ometrical approach 0. Each unitary U G S'[/(2") cor- 
responds to a point in the Riemannian manifold Af — 
SU{2^), and each Hamiltonian describes a vector in TAf, 
the tangent space of J\f. Distances on JV are defined by 
C : Af X TAf M. The minimal cost D{I, U) coincides 
with the minimal distance between / and U . 

We focus on a class of quantum control problems that 
split the space of Hamiltonians su(2") into two vector 
subspaces, I and p. Hamiltonians in [ have negligible 
cost, while those in p do not. In addition, I and p satisfy 
the set of commutation relations that define a Cartan 
decomposition (2l|: 



[[,[]C[ [p,[]=p, [P,P]C[. 



(2) 



We refer to such problems as Cartan control problems. 



Definition 1 (Cartan control problem) A Cartan 
control problem on an n-quhit system is defined as 
follows. Let [ and p be subspaces of su(2") that satisfy 
(0). Define Pi and Pp as their respective projection 
operators. The application of a Hamiltonian H for time 



dt incurs cost C{U,H) = y {H\Q H)dt, where 



g = eP+Pp, 



e< 1. 



(3) 



Given an unitary U G 5t/(2"), we wish to find H(t) on 

[0, T] that minimizes D{I, U) = /J" C{U{T), H{t))dt sub- 
ject to (QP. Alternatively, this problem can be regarded as 
computing the distance between I and U on the manifold 
N = 5J7(2") subject to the metric C. 

The 2-qubit system, where we wish synthesize U S 
SU{2") with the minimal amount of non-local interac- 
tions, is a special case of this problem. Here, I is the vec- 
tor space of single-qubit Hamiltonians and p is the vector 
space of all directions orthogonal to I. The resulting Car- 
tan control problem neglects the cost of all single-qubit 
interactions, and thus D{I, U) is a measure of the mini- 
mal amount of interactions required to synthesize U . In 
fact, it coincides with the interaction cost oiU 1^], when 
minimized over all possible drift Hamiltonians. 

The physical interpretation of n-qubit Cartan control 
problems for n > 2 is not as transparent. Although there 
exists a decomposition such that all single-qubit inter- 
actions are contained in I, [ will invariably also contain 
interactions involving an unbounded number of qubits. 
Therefore, the condition e ^ implies that in addition 
to local interactions, certain non-local interactions can 
also be applied at negligible cost. Although the solution 
for these cases does not have a direct physical applica- 
tion, other than provide a lower bound on complexity, 
it shows how geometrical methods are well adapted to 
solving a general class of problems that scale with n. 



III. SOLUTION TO THE CARTAN CONTROL 
PROBLEM 

In this section, we solve the Cartan control problem 
using geometrical methods Q- Before approaching the 
problem directly, we illustrate the intuition behind our 
approach by a simple example. 

Consider a cylindrical surface of unit radius N — 
M X [0, 2n) parameterized by standard cylindrical co- 
ordinates, z and 9 and the naturally induced metric 
CN{z,d,dz,d9) = \/dz^ + d(P. Suppose we wish to find 
the minimal distance between two points, x = (0, 0) and 
y = (0,7r/2), it is clear that geodesies between the two 
points are non-unique since the surface wraps around it- 
self. We circumvent this difficulty by introducing a sec- 
ond manifold A4 —M.'^ with the st andard Euclidean met- 
ric Cm{p, g, dp, dq) — \J dp^ + dcp- , together with a map- 
ping U : M. ^ Af oi the form [/(p, 5) — {p, g mod 27r). 
If we define [x] and [y] as the pre-image of x and y with 
respect to U. i.e: [x] = {0,2j7r}, [y] = {0, (2fc + i)7r)} 
J, fc e Z, then the distance c^Af (x, y) on Af coincides with 
the minimal distance between the sets [x] and [y] on A4 , 
i.e: tt/2. The following lemma states this more generally: 

Lemma 1 Let A4 and Af be Riemannian manifolds with 
distance measures Cm o,nd Cn. Denote the distance be- 
tween two points on A4 and Af by dM{-,-) and d]s[{-,-) 
respectively. 

Let U : A4 Af be a smooth map that preserves the 
distance, i.e: Ca/((7,v) = CN{U{q),U* (y)), where U* is 
the pushforward of U. Define [x] — {p : (t>{p) = x}, 
[y] — {q • 0(1) = y} the pre-image of x,y e A/" 
respectively. If \i £ [x] , and q G [y], then 

dhi (N, [y]) = min dM(p, q) 

qe[y] 

= min dM (p, q) = ^JV (x, y) , (4) 

pe[x\ 

where dM(p,q) denotes the distance between p and q. 
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FIG. 1: (Color onIine)If [x] and [y] in M are the pre-images 
of x and y in N, then dAf([s], [y]) = dN{x,y) provided the 
mapping U : M N is smooth. 

Proof: Suppose c?a/([2;], [y]) = k, then there exists a 
curve J C A4 that connects some p G [x] and q S [y] of 
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length k. Clearly its image r{t) = U {jit)) is a curve from 
X to y of length k in J\f. Thus (iAr(x,y) < cIm {[x], [y])- 

Now suppose dAr(x,y) — fc, such that there exists a 
curve r{t) G A/", r(0) = x, r(l) = y of length k. Given 
any p G [cc] , we show that there exists a q e [y] such that 
dN{p, q) < by constructing a curve 7 from p to q. Let 
I be a large integer and tj — jk/l, j = 0, . . . ,m, set 

7(t) = vt^. t + 7 (ij ) < i < tj+i (5) 

where m is a large integer, and vt^. satisfies 

f/X=f(*.) (6) 
Clearly in the limit n — > 00: 

Length(7(i)) = lim ^ Cm (7 (^j) , v^J = /c (7) 

Hence (iAr(x, y) > minqg^j.] (iM(p, q)- Combining the two 
results gives the desired equivalence. Symmetry implies 
dAr(x, y) > miupgjj,] dM(p, q), which establishes the de- 
sired result. ■ 

To compute distances on SU{2"), we define an Eu- 
clidean manifold Ai =M^ Denote its coordinates by 
q = {qi, . . . , (74>i_i) and tangent vectors by v. We wish 
to find a suitable metric Cm (q, v) = -y/ (vjCv) , together 
with a distance preserving map U such that Lemma [1] is 
applicable. There are many possible choices, of which we 
ideally select one where G has a simple form. 

The Cartan decomposition is one such candidate [2ll |. 
Let 3 be a maximally commuting subspace of p , then any 
unitary can be decomposed into 



(8) 



where Hi , H3 e I and H2 E The vector of matrices 
{Hi, H2, H3) completely specify U. We view this as a 
cartesian plane, (qi,q2,q3) G M, where qi is the vec- 
torization of Hi with respect to some orthonormal basis 
Bj^i, i.e. Hi = J2j and Qj^i = TT[HiBj^i]. ^ then 

defines the desired coordinate map. 

The second step is to compute the metric Cm on M 
such that CM{q,v) = CN{U{q),U* {v)). The matrix 
G can be represented by a 3 x 3 matrix of superoper- 
ators, Gij, such that a perturbation {AHi, AH2, AH3) 
has length A J^ij 
be characterized: 



V {Hi\GijHj) 



The properties of G can 



Lemma 2 Let A4 = ^ be a Riemmannian 

manifold with metric Cm{AHi, AH2, AH3) = 
AJ2y/iHi\G,jHj), and U : M ^ S'[/(2") he as 
defined by (0)! If Cniq,^) = GN{U{q),U* (v)), then G 
has the form 

eBCH^BCHL 



G 



I I . (9) 

eA{H2) + B(i/2)) 



A{H2) and B(_ff2) o,f£ H2 dependent operators that sat- 
isfy A(0) = /, B(0) = Q, {I + A(i?2)) > and 
B(i/2) > 0; and BCH denotes the Baker-Campbell- 
Hausdorff operator, which satisfies exp{— i(C + D)} = 
expl-iBCHclI?)} exp{~iD} + Oi\D\^). 



Proof: From the BCH equation, e^^e^^ = 
O(A^), thus Gij = for i ^ j. The remaining compo- 
nents can be computed by considering individual pertur- 
bations. The details are purely technical, and are left to 
the appendix. ■ 
To simplify the notation, we set L = Hi, Z = H2 and 
M ~ H3, so that a point on M is denoted by the 3- 
tuple (L, Z, M). While the explicit forms of A and B are 
complex, the metric greatly simplifies when e — > 0. 



Lemma 3 In the limit e 
unitary U is given by 



D{I,U)= min \Z 

Z:U(L,Z,M) = U 



0, the cost of synthesizing a 
(10) 



where \Z\ ~ ■\/Tr(Z^) is the trace norm of Z and U is as 
defined in (0). 







q= {L,Z,M) 
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* p = (L,0,M) 



FIG. 2: (Color online)d(p, q) = and d(0, q) 
Thus d(0,p) \Z\ as e approaches 0. 



as e ^ 0. 



Proof: Consider the triangle with vertices = (0, 0, 0), 
p = (L,0,M) and q = {L,Z,M) {Fig^. Let C be a 
constant such the operator norm |BCLl| < C and K — 
max(C, 1). The length of straight line from the origin 
to p is bounded above by eKy/ (|Lp -|- |Mp), thus so is 
d{0,p), the distance from the origin to p. Two triangle 
inequalities then bound d{0, q) from above and below: 



rf(0, q) < rf(p, q) + d(0, p) < \Z\ + e/fv/(|iP + |M|2 



d(0, q) > d(p, q) - d(0, p) > |Z| - eKy/{\L^T\MF) 

In the limit e ^ 0, d{0, q) = \Z\. Application of Lemma 
[1] gives the required result. ■ 
We now have an explicit characterization of distances 
on the coordinate manifold M. The final step is to de- 
termine the pre-image of a given imitary in S'[/(2"). We 
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use a variation of the technique developed in [l^ . Let the 
Cartan decomposition of [/ be as in ([5]), the properties 
of Cartan decompositions allow us to choose a basis such 
that the matrix representation of U can be expressed as 

U = ADB^, (11) 

where A = e'^, — e*^^ are orthogonal, and D ^ e*^ 
is diagonal [2^. is the diagonalisation of U^U and is 
hence unique up to permutation of its diagonal elements. 
We use this decomposition to find an explicit expression 
for D{I,U): 

Theorem 4 Consider the n-qubit Cartan control prob- 
lem. The minimal cost required to synthesize a unitary 
U with Cartan decomposition e'^e'^e'^^ is 



D(/,C/)= /min|eig(Z)-y|2 (12) 

V 

where C is a lattice defined by the set of points: 

£ = {(mi,TO2, . . .m2")7r : ^ = 0, Wfc G Z} (13) 

Proof: Since D is diagonal, we can describe it by a vector 
X = (xi, . . . ,X2") such that the diagonal elements of D 
take the form e"''. In addition, we know D = e'^ for 
some Z = J2j ^jBj G 3, where Bj, j = 1, ... 2" — 1 is an 
orthonormal basis for 3. 

Let A be the mapping that takes the vector x to z, the 
vector representation of Z in the Bj basis, i.e: 

eig{[AxfBk} = {xuX2,...,X2r^} (14) 

To see that A is an isometry, i.e., |.4xp — |xp, we note 
that in our particular representation (|lip . Z is diagonal, 
and thus eig { [y4x]''^Bfc} = diag { [Ax]''Bk}. 

Define Bk = diag(_B/j) as the vector formed from the 

diagonal elements of Bk, then [Ax]'' B k = x. Let B = 

[B I, B2, ■ ■ ■ , B k] be the matrix whose columns are the 
elements of Bk, then the equation can be rewritten as 
BAx = X. Since B is orthonormal, A = B~^ must be 
also, and hence preserve the Euclidean norm. So 

D(/, U) = min{|Z| : e^e^e*^ = U}, 

= min{|x| : e^ e^-^""^' e^'^ = U}. (15) 

Since permutations preserve the Euclidean norm, the 
only freedom in Xk that we need to minimize over is addi- 
tion by multiples of tt. Thus given one particular decom- 
position of a given unitary U, eig(iJ2) gives one possible 
x. The set of all vectors (permutations excluded) that 
generate U is given by {x + 1 : 1 e C}. The result fol- 
lows. ■ 
Given a unitary U E SU{2"-), the above theorem of- 
fers a systematic method to solve for the minimal cost 
required to synthesize D{I, U). 



IV. THE SINGLE QUBIT CONTROL PROBLEM 

In this section, we illustrate our result by applying 
it to the special case of single qubit optimal control. 
We wish to synthesize a particular spin i interaction 
U e SU{2). Application of magnetic fields in one par- 
ticular direction (say x) incur negligible cost, while all 
orthogonal directions require unit cost, i.e. I = Span(CT2;) 
and p = Span(tT2,j,). 

This problem is a slight variation of the single qubit 
time-optimal control problem solved in JTj. More pre- 
cisely, it corresponds to the case of a system that evolves 
under a constant magnetic field described by the Hamil- 
tonian Hd — V^cr^. We wish to synthesize a unitary U 
in minimal time, given the ability synthesize magnetic 
fields in the x direction of arbitrary strength, or reverse 
the direction of iJ^. 

Proposition 5 Let U € SU{2). Suppose we are given 
one particular decomposition 

U = exp(— ixua;) exp(— izCTz) exp(— iycTa;), (16) 

where a^, denote standard Pauli matrices, then 
D{I,U) = -^mmmez{z — 2rmr}. In particular, 

D{LU) = M forze[~n,7:]. 

Proof: Note that zaz has eigenvalues ±|, and 
C = {(mTT, -TOvr) : m e Z}. So D{I,U) = 
-\/2 miurngz \z/2 — myrp. The result follows directly. ■ 
In [ll|, a slightly different result where D{I,U) = ^ 
for z G [0, 27r] is obtained. The deviation results from 
our extra assumption that the direction of Hd can be 
reversed. The KGB result requires the unique decompo- 
sition such that 2 e [0, 27r], whereas our result applies to 
any decomposition that satisfies l|16p . 

V. CONCLUSION 

The geometrical approach provides a useful alternative 
to more algebraic methods [13, [HI, In this paper, 
we have demonstrated how we can use it to characterize 
the general Cartan control problem. In the single-qubit 
case, our result solves a slight variation of single-qubit 
time optimal control, and provides a second perspective 
to [ll| . In the two qubit case, it characterizes the minimal 
amount of non-local interactions required to synthesize a 
given interaction. 

The general n-qubit Cartan control problem that we 
have described does not have direct physical application, 
because the class of Hamiltonians assumed to be 'easy' 
to apply is too broad to be realistic. However, our re- 
sults do show an instance where the geometric formalism 
can be applied to systems of arbitrary size. By reducing 
the complex space of unitary operations into a cartesian 
coordinate system with a suitably appropriate metric. 



we circumvent much of the technical difRcuhies in alge- 
braicaUy intensive methods. 

The geometrical method outlined can convert any 
quantum control problem into a minimization of dis- 
tances between two sets in cartesian space with a suit- 
ably defined metric. This allows analytical solutions in 
special cases, such as the Cartan control problem. Al- 
ternatively there exists numerous numerical techniques 
that have been designed to solve for minimal distances 
on manifolds. Thus, the geometric formalism is a promis- 
ing method, both for solving other problems in control 
theory, and for its applications in quantum complexity. 
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APPENDIX A: DERIVATION OF THE 
COORDINATE METRIC Cm 

In this section, we provide a detailed outline of how 
the metric Cm in Lemma [2] is derived. We consider a 
point (L, Z, M) on M. and consider perturbations on each 
of L, Z and Af, which we denote by Pl, Pz and Pm 
respectively. Firstly 

U{L + APl, Z, M) = e-'(i+APL)g-^Zg-^M 

= exp[-iABCHL(PL)]C/ (Al) 
Since [L, Pl] £ I, and BCHl G [, 

{Pl\GPl) = e(PL iBCHj^BCHiPi) (A2) 
=^ Gn = eBCH^^BCHi (A3) 

Similarly, other components of G can be computed by 
perturbing Z and M: 

U{L, Z + APz,M) = e-i^e-K^+'^^'''>e''^ (A4) 
= exp [-iAe-'^Pze'^] U (A5) 
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Noting that e~^^Pze^^ £ p since [[, p] € p, we have 

{Pz\GPz) = \Pzf ^ G22=/ (A6) 
The final perturbation is slightly more complex: 
U{L, Z, M + APm) = e-~'^e"'^e"'^'^'+'^P"''> 

^ g-iLg-'*Ae-'^BCHM(PM)e'^g-iX2 
^ g-jAe-'-^e-'^BCHM(PM)e"''e"*^[/ 

_ g-i^JL,Z,M{PM)jJ (A7) 

Noting that BCHm(Pm) G P, and [p,p] C [. We can 
write: 

e-'^BCHM(PM)e^^ = a{Z)Qp{PM) + b{Z)QL{PM) 

(A8) 

Where Qp e p and Ql e [, a'^ + P = 1 and 6(0) = 0. The 
commutation relations [[, [] C I, [I, p] C p then implies 
that Jl,z,m takes the same form, i.e: 

Jl,zm{Pm) = a{Z)Qp{PM) + b{Z)QL{PM) (A9) 
Thus 

{Pm\GPm) = eaHZ) + b'iZ) G22 = I (AlO) 
In particular 

G33 |z=o=eA(Z) + B(Z) (All) 

For some positive definite operators A and B such that 
A(0) = I and B(0) = 0. Hence, in matrix representation, 
the global metric takes on a block diagonal form: 

/ eBCrf^BCHi \ 
g=\ I (A12) 

V eA{Z) + B{Z)) 
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